A new type of charged defect in amorphous chalcogenides. 
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We report on density-functional-based tight-binding (DFTB) simulations of a series of amorphous 
arsenic sulfide models. In addition to the charged coordination defects previously proposed to 
exist in chalcogenide glasses, a novel defect pair, [As4]~-[S3] + , consisting of a four-fold coordinated 
arsenic site in a seesaw configuration and a three-fold coordinated sulfur site in a planar trigonal 
configuration, was found in several models. The valence-alternation pairs [S3] + -Sj~ are converted into 
[As4]~-[S3] + pairs under HOMO-to-LUMO electronic excitation. This structural transformation is 
accompanied by a decrease in the size of the HOMO-LUMO band gap, which suggests that such 
transformations could contribute to photo-darkening in these materials. 

PACS numbers: 71.23.Cq, 61.43.Dq, 
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Amorphous chalcogenides, i.e. sulfides, selenides and 
tcllurides, are distinguished from other materials by 
their photosensitivity. These materials exhibit intriguing 
photo-induced phenomena that are not observed in their 
crystalline counterparts. Some of these unusual phenom- 
ena have many potential applications . A detailed mi- 
croscopic understanding of the origin of these phenom- 
ena 0, however, is still lacking. 

In the context of photo-induced structural changes, 
special significance is attributed to the presence of oppo- 
sitely charged coordination-defect (valence-alternation) 
pairs 3]. Normally, such defect pairs contain singly- 
coordinated chalcqgen atoms having distinct spectro- 
scopic signatures |3J. Experimentally, the concentration 
of these defects is estimated to be rather small [ij], i.e. 
10 17 cm -3 , compared with the atomic density of about 
2 x 10 25 cm~ 3 , in order quantitatively to account for the 
observed magnitude of the photo-induced effects. 

Here, we investigate the structure of, and simu- 
late photo-structural changes in, a series of computer- 
generated models of the archetypal amorphous chalco- 
genide, diarsenic trisulfide (a-As2S3). In addition to the 
well-known Q valence-alternation pairs (VAPs) [C3] + - 
Cj" and Pj-Cr, where C n /P„ stands for an n-fold coor- 
dinated chalcogen/pnictogen (e.g. sulfur/arsenic) atom, 
some of the models were found to contain the previously 
unreported [P4]~-[C3] + defect pairs. These defect pairs 
are unusual in two ways. First, there is an excess of neg- 
ative charge in the vicinity of the normally electropos- 
itive pnictogen atoms and, second, there are no under- 
coordinated atoms with dangling bonds in these local 
configurations. The latter peculiarity may be the rea- 
son why such defect pairs have not yet been identified 
exp er iment ally. 

All our simulations have been done by using a density- 
functional-based tight-binding method ||, unless speci- 
fied otherwise. We found that a basis set of s, p and d 



Slater-type orbitals for all atoms is an essential prereq- 
uisite for the observation of the defect-related effects re- 
ported here Q- Each of our models contained 60 (24 As 
and 36 S) atoms. This size of model was chosen so as 
to be big enough that the defect pairs could be accom- 
modated within the volume of the periodic simulation 
box, and yet small enough to allow for the creation and 
analysis of several statistically independent models in a 
reasonable time. The quality of the models, in terms of 
comparison with experimental neutron-diffraction data, 
is similar to that of the larger models reported in Ref. 0- 

By using NVT (constant number of particles, volume, 
and temperature) molecular-dynamics simulations, we 
first created a 24-atom model of amorphous arsenic by 
quenching from the melt. Three- fold coordination was 
imposed by iteratively modifying the nearest-neighbor 
shell of each atom and bringing the system to ther- 
mal equilibrium at room temperature. A stoichiometric 
model (model 0, in the following) of a-As2S3 was then 
created by decorating each of the As- As bonds with S 
atoms, followed by rescaling the model to the experi- 
mental density Q p = 3.186 g/cm 3 and relaxing at room 
temperature (T — 300 K). This model was then melted at 
T = 3000 K and ten more models were created by taking 
snapshots of the 3000 K molecular-dynamics trajectory 
at irregular time intervals as starting configurations, and 
following the cooling schedule described in Ref. 0. 

The defect statistics for the models containing coordi- 
nation defects are shown in Table [I] The coordination 
numbers were calculated as the number of nearest neigh- 
bors within a spherical shell of a radius corresponding 
to the position of the first minimum in the correspond- 
ing partial radial distribution function. The numbers in 
each triplet in Table [I] show, respectively, the number of 
corresponding defects in (I) the as-prepared ground-state 
optimized geometry, (II) the excited-state geometry op- 
timized under the constraint of having one electron in 
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TABLE I: Defect statistics for the models containing coordi- 
nation defects. Each triplet of numbers denotes, respectively, 
(I) the number of defects in the as-prepared ground state, (II) 
HOMO-LUMO-promoted excited state, and (III) the ground 
state resulting from the excited state. 

Number of 

Model coordination defects homopolar bonds 





As 2 


As 4 


Si 


s 3 


As- As 


S-S 





0,1,0 


0,0,0 


0,1,0 


0,0,0 


0,0,0 


0,0,0 


1 


0,1,0 


1,1,1 


0,0,0 


1,0,1 


1,1,1 


1,1,1 


2 


0,2,0 


0,0,1 


1,0,0 


1,0,1 


2,1,2 


2,2,2 


3 


0,0,0 


1,1,1 


1,1,1 


0,0,0 


1,1,1 


0,0,0 


5 


0,0,0 


0,0,0 


0,1,0 


0,0,0 


1,1,1 


1,0,1 


6 


0,0,0 


1,0,1 


0,1,0 


1,1,1 


2,2,2 


2,2,2 


10 


0,1,0 


0,0,1 


1,0,0 


1,1,1 


1,0,1 


1,1,1 



the highest occupied molecular orbital (HOMO) and one 
electron in the lowest unoccupied orbital (LUMO), and 
(III) the ground-state optimized geometry obtained from 
the excited-state configuration. Although electronic ex- 
citations where one electron is promoted from HOMO to 
LUMO Kohn-Sham states are not especially realistic, 
we simulate such excitations in order qualitatively to as- 
sess defect stability with respect to (optically-induced) 
electronic excitations. Note that in models 2 and 10, S3- 
Si defect pairs converted into AS4-S3 pairs as a result of 
the electronic excitation. 

The statistics of Mulliken charges (in atomic units) 
were calculated for all 11 models in the ground state 
before and after electronic excitation Since all extremal 
charge values, which significantly deviate from the aver- 
age values (0.32 for As, -0.21 for S), correspond to coordi- 
nation defects, these data confirm that the coordination 
defects form charged defect pairs. The maximum positive 
charge (0.46) is on a four-fold coordinated As atom that 
is part of an intimate valence-alternation pair (IVAP) 
where a singly coordinated S atom is covalently bonded 
to it. The minimum positive charge (0.07) is on an As 
atom that forms a homopolar bond with an AS4 center 
(see Table ITT)l . The minimum negative charge (-0.36) is 
on a singly coordinated S atom, and the maximum nega- 
tive charge (-0.03) is on a S atom that forms a homopolar 
bond with another S atom with a similar charge value. 

It is common Q to use the following notation for 
valence-alternation pairs: AsJ-Sj~ and S^-SJ - , where 
the superscripts show the polarity of the nominal excess 
charge on corresponding defect centers. Although the 
charges on AS4 and Si atoms are of correct polarity in an 
alloy of normally electropositive As and electronegative 
S, those on S3 are not. Therefore, it is more appropri- 
ate to use the notation [S3]" 1 " for three- fold coordinated S 
atoms, where the square brackets imply that the excess 
charge is distributed between the center and its nearest 
neighbors. Table [H] shows that there is also an excess 




(b) (c) 

FIG. 1: Planar trigonal, [Sa] + (marked by letter 'A'), and 
seesaw, [As 4 ]~ (marked by letter 'E'), configurations in (a) a 
fragment of model 2(111) (the dangling bonds show where the 
displayed configuration connects to the rest (not shown) of 
the network) and, (b) and (c), charged isolated clusters (the 
dangling bonds are terminated with hydrogen atoms). The 
shading of the As atoms is darker than that of the S atoms. 
Bond lengths are (A): (a) AB=AD=2.42, AC=2.39, EF=2.83, 
EG=2.54, EH=2.35, EI=2.33; (b) AB=2.34, AC=AD=2.33; 
(c) EF=2.55, EG=2.46, EH=2.30, EI=2.32. Color online. 



negative charge in the vicinity of AS4 centers when [S3]" 1 " 
centers are also present, and these newly-identified de- 
fects will be denoted as [As4]~ here. 

It is of interest to verify whether the excess charge has 
a similar effect on the local structure in both super-cell 
amorphous and isolated cluster models. Fig. [Ja) shows 
a fragment of model 2(111) containing both [AS4] - and 
[S3]" 1 " defects. Fig. also shows the optimized geome- 
try in (b) [SAs3H 6 ]+ and (c) [AsS4Hi]~ clusters, where 
one electron was removed from and added to the neu- 
tral system, respectively. The seesaw [As4]~ and planar 
[S3]" 1 " local configurations observed in models 1(1), 2(111) 
(see Fig.^a)), 6(1) and 10(111) are indeed very similar to 
those in the isolated clusters shown in Fig.Q(b) and (c), 
respectively. We note that similar seesaw-like structures 
are also observed in local configurations around hyperva- 
lent chalcogens, e.g. in sulfur(IV) fluoride and in certain 
tellurides llfl. and are expected from stereochemical con- 
siderations hjj . Also, the five- fold coordinated As defect 
reported in [12j can be viewed as an [AS4] ~~ center bonded 
to an additional As atom (placed on top of the structure 
depicted in Fig.^c)). Interestingly, five-membered rings 
fit conveniently to parts of the over-coordinated defect 
centers, which might contribute to the stability of this 
type of defect pair. It is also worth mentioning that 
there is a distinct pattern in the distribution of charges 
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TABLE II: Mulliken charges, in atomic units, on atoms within the [As4]~-[S3] + defect pairs. The index (I) next to the model 
number signifies the as-prepared ground-state geometry, and (III) indicates the ground-state structure obtained by geometry 
optimization of the excited-state configuration (II). The relevant fragment of model 2(111) and clusters [SAS3H6] 4 " and [ASS4H4] - 
are depicted in Fig. Q^a)-(c), respectively. 



Model 


Center 


1 




Neighb' 
2 


ors 


3 




4 


Total 


1(1) 


/A \ t~\ O -1 O 

(As) 0.313 


(S) -0.194 


(S) 


-0.281 


(S) 


-0.223 


(S) 


-0.313 


-0.698 


(S) -0.110 


/A \ 

(As) 0.363 


(As) 


i 0.379 


/A 

(As, 


i 0.362 






0.994 


2(111) 


(As) 0.329 


(As) 0.071 


(S) 


-0.221 


(S) 


-0.203 


(S) 


-0.329 


-0.353 


(S) -0.153 


(As) 0.355 


(As) 


1 0.347 


(As) 


i 0.351 






0.900 


6(1) 


(As) 0.345 


(S) -0.149 


(S) 


-0.255 


(S) 


-0.295 


(S) 


-0.225 


-0.579 


(S) -0.146 


(As) 0.375 


(As) 


1 0.376 


(As) 


1 0.384 






0.989 


10(111) 


(As) 0.326 


(As) 0.101 


(S) 


-0.242 


(S) 


-0.308 


(S) 


-0.231 


-0.354 


(S) -0.158 


(As) 0.363 


(As) 


1 0.295 


(As) 


i 0.101 






0.601 


[AsS 4 H 4 ]- 


(As) 0.280 


(S) -0.359 


(S) 


-0.320 


(S) 


-0.498 


(S) 


-0.600 


-1.497 


[SAs 3 H 6 ]+ 


(S) -0.149 


(As) 0.365 


(As) 


1 0.362 


(As) 


i 0.374 






0.952 



in the seesaw configuration, which is most clearly seen in 
the isolated [AsS-jH^] - cluster, — the more distant sulfur 
atoms F and G (see Fig.^c)) have the largest magnitude 
negative charges (see Table llTjl . 

Although the DFTB method involves a number of ap- 
proximations, it has been shown [3 that its accuracy is 
comparable to that of full DFT methods. In order to ver- 
ify whether the [As4]~-[S3] + defect pair is reproducible 
with more accurate methods, we optimized the geometry 
of model 1(1) by SIESTA O, which is a DFT method 
using an atomic-orbital basis set and pseudopotentials 
to eliminate the core electrons. Apart from insignifi- 
cant changes in bond lengths and angles, the SIESTA- 
and DFTB-optimizcd models are very similar, including 
the region with the seesaw-planar trigonal defect pair. 
We also performed high-accuracy all-electron geometry 
optimization of the isolated charged clusters shown in 
Fig. Hlb) and (c) by using NWChem HHH. The calcu- 
lations were performed at the all-electron DFT level with 
a B3LYP hybrid exchange-correlation functional and a 
6-311G** gaussian basis set which includes s, p and d 
orbitals on the sulfur and arsenic atoms. The optimized 
DFTB geometry was used as a starting configuration for 
the all-electron optimizations. After optimization, the 
shapes of the molecular clusters were preserved, apart 
from a slight increase in the bond lengths, EF=2.80 A 
and EG=2.51 A, in the "seesaw" configuration. 

We have found that the hitherto unreported [As4]~- 
defect pairs in our models are stable with respect 
to HOMO-to-LUMO electron excitations. We also per- 
formed a geometry optimization of an AS4S10H8 clus- 
ter, containing an [As4]~-[Sa] + pair by construction, 
both in the ground state and in the first singlet ex- 
cited state within the linear-response approximation to 
time-dependent density-functional theory, which gives a 
much better description of excited states compared with 



HOMO-to-LUMO electron excitations [II] (see Fig. |S). 
Although the bond AC in the ground-state structure de- 
picted in Fig. EI a) is significantly elongated, analysis of 
the electron density in the HOMO and LUMO electronic 
states (see Fig.EJb)) shows that it has a significant bond- 
ing character, implying that this cluster contains a dis- 
torted [As4]~-[S3] + pair. As seen in Fig.|3c), redistribu- 
tion of the electron density in the excited state leads to 
a symmetrization of the AS4 center and to an elongation 
of the bond BE in the region of the S3 center, where the 
LUMO electronic state is predominantly localized. Sub- 
sequent geometry optimization of the excited-state con- 
figuration results in the same ground-state geometry as 
in Fig.Efa). We observed that bond breaking/elongation 
in all our models generally occurs at the groups of atoms 
where the LUMO is localized, indicating the expected 
antibonding character of LUMO states. 

The HOMO-LUMO band-gap energies for models 0-10 
are listed in Table IIHI As in Table [IJ the three values 
correspond, respectively, to the same ground-, excited-, 
ground-state geometries. Note that in models 2 and 10, 
the magnitude of the band gap has decreased as [S3] 4 "- 
S]~ pairs were converted into [As4]~-[S3] + pairs. It should 
be noted that [S3]" 1 " centers are not necessarily conserved 
under such conversions — while in model 10 the S3 atom 
is the same in both ground states (I) and (III), in model 2 
it is not. 

In summary, we have demonstrated the existence of a 
new type of charged defect pair in amorphous arsenic sul- 
fide, namely [As4]~-[S3] + , where the As center has a "see- 
saw" configuration, and the S center is trigonal planar. A 
plausible scenario for photo-darkening in this material is 
the conversion of other types of defects containing under- 
coordinated atoms into such defect pairs due to electronic 
excitation under illumination. Interestingly, in the two 
models where such conversion has been observed, As- As 
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(b) 



(d) 



FIG. 2: An AS4S10H8 cluster containing both defect centers 
[As 4 n (marked by letter 'A') and [S 3 ] + ('B'). The shading 
of the atoms is the same as in Fig0 The black solid lines 
signify elongated bonds, (a) Optimized ground-state geome- 
try. Bond lengths are (A): AC=3.00, AD=2.32, and BE=2.4. 
(b) Isosurfaces corresponding to the value of 0.025 of elec- 
tron density in the HOMO (darker red surface) and LUMO 
(lighter cyan surface) states for the structure shown in (a) . (c) 
Optimized geometry in the first singlet excited state. Bond 
lengths are (A): AC = 2.43, AD = 2.44, and BE = 2.82. (d) 
Same as (b), but for the structure shown in (c). Color online. 



TABLE III: HOMO-LUMO band-gap energies for the all- 
heteropolar model and the models containing coordination 
defects. The indices (I), (II), (III) correspond to the triplets 
in Table □ 



Model 


(I) 


Band gap (eV) 
(II) 


(HI) 





1.46 


0.68 


1.46 


I 


1.13 


0.94 


1.13 


2 


1.41 


0.62 


1.06 


3 


1.23 


0.21 


1.24 


4 


1.47 


1.50 


1.47 


5 


1.68 


0.40 


1.67 


6 


1.47 


1.37 


1.47 


7 


1.47 


1.50 


1.47 


8 


1.38 


1.20 


1.38 


9 


1.65 


1.31 


1.65 


10 


1.20 


0.85 


1.02 



homopolar bonds were formed as part of the [AS4] ~ cen- 
ters. Given that their bond lengths are greater than the 
average bond length in the material, such bond formation 
could also contribute to photo-expansion of the material. 
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